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ABSTRACT 

We investigate bulge and disk scaling relations using a volume-corrected sample of early- to 
intermediate-type disk galaxies in which, importantly, the biasing flux from additional nuclear com- 
ponents has been modeled and removed. Structural parameters are obtained from a seeing-convolved, 
bulge+disk+nuclear-component decomposition applied to near-infrared surface brightness profiles 
spanning ~10 pc to the outer disk. Bulge and disk parameters, and bulge-to-disk ratios, are ana- 
lyzed as a function of bulge luminosity, disk luminosity, galaxy central velocity dispersion, and galaxy 
Hubble type. Mathematical expressions are given for the stronger relations, which can be used to 
test and constrain galaxy formation models. Photometric parameters of both bulges and disks are 
observed to correlate with bulge luminosity and with central velocity dispersion. In contrast, for the 
unbarred, early to intermediate types covered by the sample, Hubble type does not correlate with 
bulge and disk components, nor their various ratios. In this sense, the early-to-intermediate spiral 
Hubble sequence is scale-free. However, galaxies themselves are not scale-free, the critical scale be- 
ing the luminosity of the bulge. Bulge luminosity is shown to affect the disk parameters, such that 
central surface brightness becomes fainter, and scale-length bigger, with bulge luminosity. The lack 
of significant correlations between bulge pararmeters (size, luminosity or density) on disk luminosity, 
remains a challenge for secular evolution models of bulge growth. 
Subject headings: galaxies : spiral — galaxies : structure 



1. INTRODUCTION 

In a companion paper (Balcells, Graham, & Peletier 
2007, hereafter Paper III), we present high-resolution 
near-infrared (NIR) surface brightness profiles of carly- 
to intermediate-type disk galaxies. These are derived 
from Hubble Space Telescope (HST) NICMOS images 
and extended with NIR ground-based profiles from 
UKIRT images. That paper provides profile decomposi- 
tions into bulge, disk, and nuclear components, and ana- 
lyzes the nuclear properties of the bulges. In the present 
paper, we use the profile decompositions of Paper III to 
analyze the global scaling relations for bulges and disks, 
free from the biasing influence of additional nuclear com- 
ponents. 

Global scaling relations provide useful diagnostics on 
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the structure of disk galaxies. For bulge components, 
comparison with the scaling relations of elliptical galaxies 
shows to what degree bulges are either similar or differ- 
ent to ellipticals. Whether bulges resemble spheroids or 
disks is still an open question (Wyse, Gilmore, & Franx 
1997; Kormendy & Kennicutt 2004). To give one exam- 
ple of the complexity of this issue, some apparent bulges 
in face-on, barred SO galaxies are in fact inner disks (Ko- 
rmendy 1993; Erwin et al. 2003), which of course has im- 
plications for galaxy formation models. For both bulges 
and the outer, large-scale disks, the scaling relations may 
provide clues on the long-standing question of the origin 
of the Hubble sequence. 

The sample of de Jong & van der Kruit (1994) has 
provided a useful reference for intermediate- to late-type 
field disk galaxies, and its scaling relations have been 
analyzed by de Jong (1996) and Courteau, de Jong, 
& Broeils (1996). Graham (2001, hereafter G01) and 
MacArthur, Courteau, & Holtzman (2003) provide a 
more recent study of intermediate- to late-type spirals, 
using Sersic fits for the bulge. Hunt, Pierini, & Giova- 
nardi (2003, hereafter HPG03) have also analyzed bulge 
and disk scaling relations from a sample of disk galax- 
ies in the Perseus-Pisces supercluster; their sample com- 
prises mostly Sb and later- type spirals. At the other 
end of the Hubble sequence, SO galaxies are often dis- 
cussed together with ellipticals, but are rarely compared 
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to spiral galaxies. The study presented in this paper is 
complementary to the ones above as it samples early- to 
intermediate-type, SO through Sbc, disk galaxies. More- 
over, the analysis presented here uses HST images to 
subtract the flux of additional nuclear components such 
as star clusters and nuclear disks, which have biased pre- 
vious studies of bulge, and hence disk, parameters. This 
is perhaps best evidenced through the reporting of 14- 
like bulges with Sersic indices n ~ 4 — 6 when using 
low-resolution ground-based data which smear out the 
flux from the unresolved nuclear components (e.g., An- 
dredakis, Peletier, & Balcells 1995, hereafter APB95). 
Higher-resolution studies with HST have since revealed 
that the majority of such galaxies have noticeably less 
concentrated bulges, with n < 3, but clear additional 
components (Balcells et al. 2003, hereafter Paper II; Pa- 
per III). 

Bulge and disk scaling relations have traditionally 
been studied as a function of Hubble type, commonly 
parametrized with the revised type index T (de Vau- 
couleurs et al. 1991, hereafter RC3). The lack of a cor- 
relation between the ratio of the bulge effective radius 
R c and the disk scale length h with Hubble type, led to 
the statement that the spiral Hubble sequence is scale- 
free (de Jong 1996; Courteau et al. 1996). In a scale- 
free situation, an 'ice-berg' model, in which bulge sur- 
face brightness rather than effective radius determines 
how much the bulge protrudes above the disk, might ex- 
plain the higher prominence of earlier-type bulges (G01). 
That the Hubble sequence is scale-free has implications 
for bulge formation models, as secular evolution models 
(Pfenniger & Norman 1990; Norman, Sellwood, & Hasan 
1996) predict that bulge and disk scale lengths are cor- 
related (Combes et al. 1990). 

The parameters used for the scaling relations, derived 
from NIR images, are less affected by dust extinction 
than those derived from optical images. Furthermore, 
M I L variations with population age and metallicity are 
small in the NIR, yielding small differences between the 
photometric length-scales and the stellar-mass scales of 
the galaxies. We will show that bulge luminosity is in- 
deed a key yardstick that traces the values of bulge, disk, 
and global galaxy parameters. Section [2] describes the 
sample and the profile decomposition. Section 12.11 dis- 
cusses selection biases and provides a volume correction. 
The main results of the paper are contained in Figures [TJ 
[2j [3l and [4] and analyzed in Section [3l for bulges ( 33. ip . 
disks ( 33.2j) . and bulge-to-disk scaling relations ( §3.3|) . 
We discuss two specific issues in Section HJ the trends 
of the galaxy central surface brightness with spheroid 
luminosity (§ I4.1|) . and the scale-free nature of the Hub- 
ble sequence (§ 14. 2[) . A Hubble constant of Hq = 75 
kms _1 Mpc _1 is used throughout. 

2. GALAXY SAMPLE AND DATA 

The galaxy sample is described in detail in Paper III. It 
comprises 19 galaxies of types SO-Sbc extracted from the 
diameter-limited sample of inclined galaxies of Balcells 
& Peletier (1994). The diameter limit puts the selected 
galaxy diameters in the broad range 15 kpc to 80 kpc. 
The selection excluded very dusty bulges, so the sample 
has a slight bias toward quiescent bulges. Basic prop- 
erties of the sample are listed in Table 1 of Paper III. 
The 19 galaxies were observed with NICMOS on HST 



through the F160W filter. Data reduction is described in 
Peletier et al. (1999, hereafter Paper I). The derivation of 
surface brightness profiles is described in Paper III. K- 
band surface brightness profiles obtained from UKIRT 
images (APB95; Peletier & Balcells 1997) were scaled 
to the H-band and linked smoothly to the HST profiles, 
yielding surface brightness profiles which span from ~20 
pc to several kpc, thus covering the nucleus, bulge, and 
disk-dominated region of each galaxy. 

We performed a one-dimensional profile decomposition 
using an exponential model for the disk and Sersic (1963) 
model (Graham & Driver 2005) for the bulge, 

I(R) = 1(0) exp{-6„ (R/R c ) 1/n }, (1) 

where R c encloses half the model light, n measures the 
curvature of the profile, and b n w 1.9992n — 0.3271. 
The decomposition is described in detail in Papers II 
and III. Due to the presence of positive nuclear residuals 
in some galaxies, unresolved (PS) or resolved (exponen- 
tial) components were added to the fitting function. The 
properties of these nuclear components are analyzed in 
Paper III, which also shows the profiles, the fits, and 
lists the best-fit parameters (Tables 2-4 of that paper) . 
Those parameters provide the basis for the analysis of 
the scaling relations presented here. Our analysis ignores 
galaxian sub-components that may be found outside the 
bulges, such as rings and lenses (Prieto et al. 2001). We 
also ignore bar components. Our sample selection took 
galaxies listed as unbarred, although the presence of bars 
in some of the more edge-on cases cannot be ruled out. 

2.1. Sample completeness 

Given the small sample size (N = 19) and the various 
selection processes involved, we have checked to what 
degree our sample is a fair representation of the local SO- 
Sbc galaxy population. We follow the standard V/V meix 
formalism developed by Thuan & Seitzer (1977). Briefly, 
for a diameter-limited sample, under the assumption that 
surface brightness is independent of distance, a galaxy's 
angular diameter 9 is inversely proportional to its dis- 
tance, and it is straightforward to compute the maximum 
distance d max at which the galaxy would still be included 
in a sample limited by 9 > 0l, with 9^ — 2' in our case. 
We then compute the volume V of a sphere of radius 
equal to the distance d to the galaxy, and the volume 
Knax of the sphere out to d max . The ratio of volumes is 
V/V max = (0l/0) 3 - For objects randomly distributed in 
space, V/V max should be uniformly distributed between 

and 1, with a mean (V/V max ) = 0.5 ±l/y/[l2 N), where 
N is the number of objects. 

For the original sample selected in Balcells & Peletier 
(1994), comprising 43 objects, the UGC red diameters 
yield (V/V max ) = 0.480 ± 0.044, hence that sample is 
statistically complete. The sample of 19 objects imaged 
with HST has (V/V max ) = 0.413 ± 0.066. This shows 
that the HST sample is mildly biased toward nearby ob- 
jects; this reflects the selection of targets for HST imag- 
ing, which, other things being equal, favored high spatial 
resolution. The bias introduced is small, in any case. 

By weighting each galaxy by 1/V^ lax , our diameter- 
limited sample mimics a volume-limited sample. The 
weighting is essential whenever mean values or volume- 
related quantities are sought; we have also explored its 
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use when computing scaling relations between various 
galaxy parameters, see $3] V max is given by 

4tt n3 _ 4tt f d x 6 \ 3 

''max — ""jj~("max) — I ^ I • (,2) 

The distribution of 1/V max with bulge if -band abso- 
lute magnitude will be shown in Figure[T]. 

3. GLOBAL PARAMETERS 

This section presents global scaling relations inferred 
from the bulge-disk decompositions. Galaxy param- 
eters are corrected for foreground Galactic extinction 
(Schlcgcl, Finkbciner, & Davis 1998), (l + z) 4 cosmolog- 
ical dimming, and iVcorrection (Poggianti 1997). The 
disk parameters have been corrected to a face-on aspect, 
assuming transparent disks. The corrected parameters 
are listed in Table 3 of Paper III. FiguresHHH display 
the dependencies of global parameters on bulge absolute 
magnitude, disk absolute magnitude, central velocity dis- 
persion, and galaxy Hubble type, respectively. Strong 
correlations with bulge absolute magnitude are generally 
found, with one single deviant point at the faint end. 
This corresponds to NGC 5577 (Sbc), and we exclude 
this galaxy from the regression analysis in order not to 
bias the global trends. We computed the scaling relations 
via unweighted orthogonal regressions, following the al- 
gorithm of York (1966). We explored weighted orthog- 
onal regressions using 1/Vmax weights: resulting fitting 
coefficients changed by typically ~10%, but a few of the 
weighted regressions failed to converge, hence, for the 
sake of a uniform treatment, the parameters presented 
in the following subsections are those of the unweighted 
regressions. For each correlation we list the Spearman 
rank-order correlation coefficient (Sr) and the probabil- 
ity of a null correlation (P nu n). 

3.1. Bulge parameters 
3.1.1. Sizes 

Bulge effective radii (Fig.[T^,) range from 0.1 to 1 kpc. 
They increase with bulge luminosity, as found for ellip- 
ticals (e.g. Hubble 1926; Binggeli, Sandage, & Tarenghi 
1984, hereafter VCC). An orthogonal regression to the 
R e ~LK,Bui relation gives 5 

i? c /kpc = 1O- o - 52±o - O4 (Lk,Bui/1O 1o Lk, ) o - 39±0 - 06 (3) 

(Sr = —0.87; Pnuii = 3.0- 10~ 4 ). A strong correlation ex- 
ists also between R c and bulge central velocity dispersion 
(Fig.Eh.; Sr = 0.74; P nuU = 2.3 • 10" 3 ). In contrast, R c 
shows no dependence on T (Fig. [4^; Sr = —0.17; P nu n 
= 0.46). Both of these results are similar to those of 
HPG03. The correlation of R D with disk absolute mag- 
nitude is also null (Fig.UJi; Sr = -0.32; P null = 0.18). 

Mollenhoff & Heidt (2001) report R c ~ i 84 for their 
bulge sample. Part of the difference between the ex- 
ponents appears to be due to the type of regression per- 
formed, as an orthogonal regression to their K-band data 
yields R c ~ L°- 51 ± o oe ^ j. e ., at 2-er from our result. Our 
two samples cover similar bulge luminosity ranges, but 
their R c go up to 10 kpc for the largest bulges. Their 

5 M KjQ = 3.41 (Allen 1973) 



fits, using ground-based data, include any nuclear com- 
ponents as part of the bulge. Given the 90% detection 
frequency of nuclear components in our sample, such 
components are likely to be present in their galaxies, 
which may boost Sersic indices up to n ~ 6 (Paper II); 
and such biases yield higher values of R e . 

For elliptical galaxies, VCC show that the R c ~ L a 
relation has a < 0.5 for low-luminosity ellipticals (Mb + 
5 log/150 > —20.5), and a > 0.5 for giant ellipticals 
(Mb + 5 \ogh 50 < —20.5). Our bulges have absolute blue 
magnitudes in the range —14 < Mb + 5 log h 50 < —20 6 . 
Over that range, the R e — L relation for VCC ellipticals 
is approximately R c ~ L 0A5 . Comparing with eqn.[3l we 
conclude that the size-luminosity relations for bulges has 
similar slope to that of ellipticals of the same luminosity. 

3.1.2. Effective surface brightness 

The mean effective surface brightness is readily ob- 
tained from R c as (a) e = K + 5 log R c + 2.5 log 2ir, 
where K is the total apparent iVband magnitude of the 
bulge. VCC showed that, for faint ellipticals, (fi) e be- 
comes brighter with luminosity, while it becomes fainter 
with increasing luminosity for giants. For bulges, the 
R c — L relation (eqn.[3]) indicates that (/i) e brightens 
with bulge luminosity over the bulges' luminosity range. 
The R c ~ L 5 which corresponds to (/i) c being indepen- 
dent of luminosity isshown with an offset dashed line in 

Fig.no. 

The effective surface brightness, fi c = fi(R c ) = 
—2.5 log I c , is shown against TVband absolute bulge 
magnitude in Figure[T}D. An orthogonal regression gives 

W-&x,o = lO- 4 - 78±o - 11 (Ljr,Bui/lO lo ijc,0)°- 3a±0 - 19 . 

(4) 

(Sr = 0.34; P nu ii = 0.16). Here and in following sub- 
sections, if-band surface brightness are derived from H- 
band values using H — K — 0.23. The distribution is 
nearly flat, as are the trends of \x e with the disk lumi- 
nosity Lif.oisk (Fig.[2b; S R = 0.41; P nu n = 0.087), with 
log(oo) (Fig.Eb; Sr = -0.48; P nuU = 0.065), and with T 
(Fig.Hb; Sr = 0.43; P nul] = 0.069). 

3.1.3. Central surface brightness 

Of key interest are the trends of central surface bright- 
ness with luminosity. If bulge surface brightness pro- 
files were homologous, i.e., if they had the same profile 
shapes, the trends found in §3. 1.21 for (/i) c and /x e would 
imply central densities nearly independent of luminosity. 
However, bulge profiles are not homologous, and, hence, 
peak surface brightness needs to be measured closer to 
the centers (Binggeli & Cameron 1991). 

Even with the high spatial resolution of the HST, the 
various approaches that may be envisaged to define and 
measure central surface brightness lead to different re- 
sults. Each approach has its own merits. One measure 
is provided by the inwardly-extrapolated central value 
of the bulge Sersic profile, /Uo,Ser- Using an extrapo- 
lated central value is common practice for exponential 
and other profile models that have finite central density, 
and so we use this approach for our Sersic spheroids. As 

6 Given a mean color of B — K = 4.0 for our bulges (Peletier & 
Balcells 1997). 
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Fig. 1. — The dependence of bulge and disk parameters, and bulge-disk ratios, on the A"-band bulge absolute magnitude. All parameters 
are from Table 3 in Paper III, except for the velocity dispersions which are from Falcon-Barroso et al. (2002), and are listed in Table 1 
in Paper III. Solid lines show orthogonal regressions to the data points; the corresponding relations are given in |J3] The open circle 
corresponds to NGC 5577, an outliner in most of the distributions; it has been excluded when computing all of the regressions, (a) Effective 
radius of the (Sersic) bulge component. The dashed line shows L ~ R c 2 , offset for clarity, (b) Effective surface brightness of the Sersic 
component, (c) Extrapolated _ff-band central surface brightness of the bulge component, (d) Disk major-axis scale length, (e) Face-on 
H-band extrapolated disk central surface brightness. The horizontal dotted line is the canonical Freeman value, using the mean color for 
the sample B — H = 3.7. (/) Central velocity dispersion. The dashed line gives the Faber-Jackson relation for Coma ellipticals, from Pahre 
et al. (1998). (g) Bulge-to-disk central brightness ratio log(7o b/Io d)- (h) Ratio R c /h between the bulge effective radius and the disk 
scale length, (i) Bulge-to-disk luminosity ratio B/D. (j) 1/Vmax volume correction, normalized to the maximum 1/V^nax- 
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Fig. 1. — See caption in previous page. 



compared to a direct reading of the fi(r) profile, no,Ser 
avoids the biases introduced by additional components 
present in the galaxy nuclei (Paper III), and hence may 
be seen as a more accurate estimate of the central sur- 
face brightness of the bulge component. We will show 



here that ^o,Ser shows a strong, monotonic trend with 
luminosity for bulges, as previously known for dwarf and 
giant ellipticals. 

Figure[TJ; shows /io,Scr = —2.5 log io,Scr against the 
bulge absolute magnitude. The distribution shows that 
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Fig. 2. — Identical to Figure[T] except for the use of the _K"-band disk absolute magnitude in the abscissae. 



the central surface brightness of bulges is a steep function 
of luminosity (compare to Fig.[T})). After offsetting /Uo.Ser 
(iJ-band) to the if-band with /^o,Ser,A' = A*o,Scr — 0.23, 
an orthogonal regression to the data points gives 

T IT nn -3.37±0.15/ r /iniOr \0.81±0.33 

(5) 

(Sr = 0.49; P nu ii = 0.044). The relation is close 
to the maximum slope that is compatible with an in- 
crease of R e with spheroid luminosity (given Lsersic = 
1(0) R c 2 2irnr(2n)/b% a ). There is no sign of a turn-over 
at high luminosities. Indeed, if n increases with lumi- 
nosity, then the effective surface brightness turn-over at 
bright magnitudes does not imply a turn-over of cen- 
tral surface brightness (Binggeli & Jerjen 1998; Graham 
&; Guzman 2003). Rather, fig continues to rise with 



luminosity due to the higher n from higher luminosity 
spheroids. The dependence of /zo,s e r on T is not signif- 
icant (Figure|4j;; P nu u = 0.97), indicating that Hubble 
type index does not determine central galaxy density. 
Similarly, /io,Scr does not correlate with the disk lumi- 
nosity (Figure^; P nu u =0.57). 

In Figured we plot the Sersic-extrapolated central 
surface brightness /io,Scr of our bulges against absolute 
magnitude, together with ellipticals from Caon et al. 
(1993) and dwarf ellipticals from Binggeli & Jerjen (1998) 
and Graham & Guzman (2003). The values of /io,Scr from 
bulges follows a general trend of brighter ^o,Sor for more 
luminous objects, common to spheroids of vastly different 
luminosities: objects plotted in Figure[5^i include dwarf 
ellipticals and giant, 'core' ellipticals present in the Caon 
et al. sample. (At the faint end, bulges are brighter than 
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Fig. 3. — Identical to FigureJT] except for the use of the central velocity dispersion in the abscissae. 



dEs of the same luminosity; we argued in Paper III that 
those nuclear disk components that do are not identified 
by our profile decomposition, might contribute to this 
offset.) 

The quantity /io,sor is a useful parametrization of a 
bulge's central surface brightness, but it can differ from 
the galaxy's central surface brightness, to which the nu- 
clear components also contribute. The latter quantity is 
relevant as a diagnostic of dissipation during galaxy for- 
mation, and, e.g., for inferences of the nuclear evolution 
during galaxy mergers. Because most of the nuclear com- 
ponents are unresolved even by HST (Paper III), the true 
galaxy's central surface brightness is generally not acces- 
sible from the present data. We have therefore explored 
the trends defined by various measurements of n(r) at a 
fixed, small angular radius; linear radius; or a fraction 



of the bulge's R c . The various measurements yield the 
same trends, namely, a brightening with increasing bulge 
luminosity. As an example, Figure[5j3 shows the surface 
brightness measured at a radius of 50 pc against absolute 
magnitude for our bulges, together with those of ellipti- 
cals (Lauer et al. 1995, hereafter L95; Caon et al. 1993), 
and dwarf ellipticals (Binggeli & Jerjen 1998, Virgo clus- 
ter data; Graham & Guzman 2003, Coma cluster) 7 . The 
choice of radius is driven by the angular resolution of 
WFPC2 and the distance to the Coma cluster. The dis- 

7 Surface brightness values from the literature have been scaled 
to the //-band as follows: for ellipti cals, we use V — H values 
derived from NED (http://nedwww.ipac.caltech.edu/ I when avail- 
able, otherwise we set V — H = 3.U. For dwarf ellipticals, we use 
V — H = 2.5, which corresponds to a stellar population of metal- 
licity 0.4 times Solar, age 5 to 10 Gyr and Salpeter IMF, from 
Vazdekis et al. (1996). 
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tribution in Figured is similar to those of e.g. Phillips 
et al. (1996, their Figure 6); Faber et al. (1997, hereafter 
F97, their Figure 4); Graham & Guzman (2003, then- 
Figure 9). Both for bulges and 'power-law' ellipticals, 
brightness increases monotonically with luminosity, the 
trend only breaking for the 'core' ellipticals due to the 
flattening of their profiles in the inner 100-200 pc (L95; 
Graham et al. 2003; Trujillo et al. 2004). We find a sim- 
ilar trend of brighter surface brightness with luminosity 
for any measurement of /Lt cen down to the resolution limit 
of our data, 10-20 pc. 

3.1.4. Bulge concentration 

Figure[6] displays the behavior of the bulge Sersic in- 
dex n, a measure of how centrally concentrated the stel- 
lar distribution is (Trujillo et al. 2001b) — as a func- 



tion of galaxy Type, log(B/D), and Mk,Bu\.- We find 
values of n from 0.7 to 3.1, a lower range than typi- 
cally reported from fits to bulge surface brightness pro- 
files of early-type disk galaxies using ground-based data 
(APB95; Khoshroshahi et al. 2000; Mollenhoff & Heidt 
2001). Paper II discusses that the higher values of n de- 
rived from ground-based data come from the inclusion 
of light from additional, distinct, central components, 
which forces the Sersic index up. Once the central com- 
ponents are dealt with, the bulge light distribution shows 
values of n lower than the classical de Vaucouleurs 14 
behavior. The value of n increases toward earlier types 
and toward higher values of B/D (Fig.[S^,b), as found 
by previous works (APB95; G01; Trujillo et al. 2001a). 

Virgo cluster ellipticals from Caon et al. (1993), and 
Coma cluster dwarf ellipticals from the -ffST-based study 
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Fig. 5. — The central //-band surface brightness of spheroids 
vs. their B-band absolute magnitude, (a) Extrapolated central 
surface brightness of the Sersic fit to the profile. (6) Measured sur- 
face brightness at 50 pc from the center. All values are corrected 
for Galactic extinction, cosmological dimming and K-correction. 
Filled circles: bulges, this work. Asterisks: ellipticals from Caon 
et al. (1993). Triangles: "power-law" ellipticals from L95. Open 
squares: "core" ellipticals from L95. Open stars: Coma dwarf el- 
lipticals from Graham & Guzman (2003). Points: Virgo dwarf 
ellipticals from Binggeli & Jerjen (1998). See footnote 7 for con- 
versions. 



of Graham & Guzman (2003) are included in Figure^ 8 , 
as well as late-type bulges from de Jong (1996) and ana- 
lyzed by G01. Our bulges follow the sequence defined by 
the Caon et al. ellipticals, while the later-type bulges of 
G01 are offset to lower values of n. The reason for this 
offset is unclear; it could be related to the differences 
in viewing angles between the face-on de Jong sample 
and our inclined sample (see also Mollenhoff et al. 2006). 
Also, it is plausible that our Sersic fits of the four lowest 
luminosity, inclined bulges are affected by undetected in- 
ner exponential components — raising the bulge shape 
index above its true value. 

The dwarf ellipticals also have overall higher values of 
n than the late-type bulges of similar magnitude. This 
difference appears also for the dE sample of Binggeli & 
Jerjen (1998) and is probably real. The offset is not 

8 A constant color term _F606W — K = 2.7 was used to scale the 
Graham & Guzman dE absolute magnitudes, which corresponds 
to a stellar population of metallicity 0.4 times Solar, age 5 to 10 
Gyr and Salpeter IMF, from Vazdekis et al. (1996). 



due to an incorrect color correction for the dEs; a lower 
metallicity would move the dE total luminosities faint- 
ward, while solar mctallicitics would make them brighter 
by only ~0.3 mag. It is unlikely that nuclear sources in 
the dEs are biasing the Sersic fits, as Graham & Guzman 
allowed for additional nuclear components in their fits, 
as we have done with the bulges sample. 

3.1.5. The Faber- Jackson relation for Bulges 

We show the Faber & Jackson (1976) L — a$ re- 
lation (see also Poveda 1961; Minkowski 1962; Fish 
1964) for bulges in FigureQJ. The velocity dispersions 
are aperture-corrected values from Call triplet spectra 
(Falcon-Barroso et al. 2002). We find a tight correlation 
which yields 



a = 10 



2.04±0.04 



-0.84; P„ 



(L K , Bul /10 w L K , e y-^™, (6) 
= 5.6 ■ 10 -4 ) giving a relation of 



(Sr 

L ~ <7q' 9±0 ' 5 . Our exponent is lower than that found for 
luminous ellipticals: Faber & Jackson give L ~ a , while 
Pahre et al. (1998) give L K ~ cr 41± ° r for Coma ellipti- 
cals. The latter relation is plotted in Figure[Tf, (dashed 
line). Formally the two relations differ. However, the 
slope for Pahre et al.'s sample is largely determined by 
objects in the luminous range Mk < —22 mag, where 
few bulges lie. The deviations onset at Mk > — 22 mag, 
where the Pahre et al.'s data show much scatter (sec 
Fig. 2 of Falcon-Barroso et al. 2002). In our domain, the 
slope defined by bulges is very well constrained despite 
our small sample size. It is plausible that, if Pahre et al.'s 
dispersions come from major axis spectra, rotation might 
lead to overestimating the velocity dispersions. We note 
that the FP of our bulges sample is very close to the 
Jorgensen et al. (1996) FP of Coma ellipticals (Falcon- 
Barroso et al. 2002). 

The shallower slope we find is in agreement with the 
studies of faint elliptical galaxies with comparable lumi- 
nosities. Tonry (1981) repored a faint-end slope of ~3, 
while Davies et al. (1983) and Held et al. (1992) report a 
slope of ^2.5 for ellipticals fainter than -20 B-magnitudes 
and dwarf ellipticals, respectively. The recent analysis of 
Coma faint early-type galaxies by Matkovic & Guzman 
(2005) gives a value of 2.01 ± 0.36. 

Interestingly, central velocity dispersion shows no de- 
pendency with galaxy T (Fig.|4f; Sr = —0.38; P nu n = 
0.11), indicating that the bulge central potential does not 
vary along the early-to-intermediate Hubble sequence, 
except, perhaps, for our latest (Sbc) galaxies. Hence, 
along this part of the Hubble sequence, Hubble type is 
not determined by the depth of the central potential of 
the galaxy. 

3.2. Disk parameters 

Our analysis of disk parameters differs from similar 
studies in the literature, in that we additionally focus 
on the scaling of disk parameters with bulge luminos- 
ity (Figs.[TH,e). Such relations may provide useful tests 
of hierarchical galaxy formation models (e.g. Cole et al. 
2000; Navarro & Steinmetz 2000), given that, in those 
models, disks largely grow around pre-existing bulges of 
merger origin. To facilitate comparison with other works, 
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Fig. 6. — The bulge Sersic index n is plotted against: (a) the revised morphological type index T from the RC3; (6) B / Dderived from the 
best-fit parameters; and (c) the bulge i\-band absolute magnitude derived from B/Dand the galaxy i^-band absolute magnitude, corrected 
for Galactic extinction, cosmological dimming and K-correction. Filled circles: bulges, this work. Crosses: bulges from the de Jong & van 
der Kruit (1994) sample, as analyzed by G01. Triangles: Coma dwarf ellipticals from Graham & Guzman (2003). Squares: Virgo ellipticals 
from Caon et al. (1993). 



we also show the trends of disk parameters with disk ab- 
solute magnitude (Figs. EH, e), with central velocity dis- 
persion (Figs.[5Ji,e), and on Hubble type T (Figs.|4ji,e). 
Most disk parameters scale with bulge luminosity, and 
do not with Hubble type, as noted previously by oth- 
ers. Perhaps more interestingly, we will show that, for 
the types analyzed in the present study, disk parameters 
show weak or null correlations with K-band disk lumi- 
nosity. 

The disk scale length h increases with bulge luminosity 
(Fig.rjji) as 

tykpc = 1O o - 35±o - o4 (L a -,bui/1O 1o L/<, )°- 36±0 - 08 . (7) 

(Sr. = 0.84; P nu n = 5.1 • 10 -4 ). It has a similar depen- 
dency on bulge luminosity than the bulge R c does, which 
leads to the result that the ratio R /h must be indepen- 
dent of bulge luminosity, which we discuss in §3.31 

As expected, h increases with disk luminosity Lk,Dis 
(Fig.[2]i). However, the relation, 

Vkpc=lO OO34±olu (iK,Di 8k /lO 11 ^,0) O55±014 . (8) 

is weaker (5r = 0.62; P nu ii = 1-1 • 10~ 2 ) than found 
against L^ B ul- 

The dependence of h with Hubble type index T is 
shown in Figure^!. Disk scale lengths show no trends 
with T (S R = 0.16; P nuU = 0.49), a fact already noted 
by Freeman (1970). This indicates that the increase of 
h with bulge luminosity (Fig.QJl; eqn.[7]) is not a conse- 
quence of earlier-type galaxies in our sample being intrin- 
sically larger galaxies. Disk scale lengths are also rather 
constant with Hubble type in the samples of HPG03 (H- 
band data) and of de Jong (1996, if-band), which are 
overall later type than ours. 

The disk face-on central surface brightness /uo,Disk = 
—2.5 log /o.Disk (Fig.QJ) also scales with bulge luminos- 
ity: 



T IT 1 n-5.55i0.06/ r /inlOr \-0.28±0.( 

i.e., a gentle but well-defined faintward trend with in- 
creasing bulge luminosity (Sr = —0.60; P nu n = 0.013). 
For reference, the Freeman (1970) canonical disk central 
surface brightness is shown (horizontal dashed line) . The 
disk fiQ does not vary with T over our range of Hubble 
types (Fig Us; S R = 0.13; P nuU = 0.57). G01 and HPG03 
find a similar behavior, together with a mild decrease for 
later types (T>4). 

The disk central surface brightness is shown against 
disk absolute magnitude in Figure^. We see no correla- 
tion, and one may be tempted to see this distribution as a 
manifestation of Freeman's law; in fact, comparison with 
Figure[T] (^o.Disk against Mk,Bui) shows that a second 
parameter, namely Mx,Bub is responsible for the width 
of the distribution, and that the lack of a correlation 
is solely a consequence of a choice of independent vari- 
able, Af^.Disk, that scrambles the distribution of /Uo,Disk 
against Mrt.buI- 

3.3. Scaling of bulge and disk parameters 

Figures[Tg-i,[3g-i, and[4g-i show how the bulge-to-disk 
ratios of central brightness, spatial scales and luminosi- 
ties depend on the bulge absolute magnitude, central ve- 
locity dispersion, and Hubble type index T, respectively. 
The dramatic rise of log(/o,Bui/^o,Disk) with bulge lumi- 
nosity (10 1 - 10 3 , Fig.QJ; Sr. = -0.65; P nuU = 7.6 • 10~ 3 ) 
is dominated by the brightening of the bulge central in- 
tensity with increasing bulge Sersic index n. In contrast, 
the position along the Hubble sequence has a weak rela- 
tion with log(/o,Bui/^o,Disk): figureSg shows a big scatter 
with a drop at only the latest (Sbc) type (Sr = —0.13; 
PnuU = 0.57). A null trend is found also with disk abso- 
lute magnitude (Fig.[2g; P nuU = 0.42). 

The ratio of spatial scales R e /h is shown against bulge 
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absolute magnitude, central velocity dispersion and T in 
FiguresQJi, [3ji, Eh, respectively. R c /h shows no definite 
trend with any of these parameters: P nu \\ = 0.41, 0.76, 
and 0.48 for Mk,Bu1, log(cro) and T, respectively. Similar 
dependencies of R c /h with n and with T have been found 
for later-type spirals (de Jong 1996; G01; HPG03). The 
non-dependence of R c /h with T is the basis for the state- 
ment by Courteau et al. (1996) that the Hubble sequence 
is scale-free, a topic we address in §4.21 

Finally, as expected, B/D increases with bulge- 
luminosity (Fig.QJ). The relation is 

B/D = lO-°- 59±olo (i A ',Bui/lO lo iK, ) o - 62±015 , (10) 

(Sr = -0.70; P nuU = 3.7- 10" 3 ). The increase in B/D is 
dominated by the increase in bulge luminosity, as already 
pointed out by Trujillo et al. (2002). B/D increases by a 
factor of nearly 100 over 8 mag in Mk,Bm\- In contrast, a 
null trend is found with disk absolute magnitude (Fig.[2j; 
P nuU =0.86). 

Equation 1101 implies that disk luminosities increase as 
^Disk ~ -^Buige ' 38 - B/D also increases with log(<7o) 
(Fig.H; S R = 0.69; P nun = 6.2 -IO" 3 ), and B/D depends 
only weakly on T (Fig.ffl; S K = -0.48; P nu n = 0.041, 
driven by the drop at T = 4), and with considerable 
scatter (see also de Jong 1996), highlighting the problem 
of using T-based B/D ratios (e.g., Baggett, Baggett, & 
Anderson 1998). Although B/D appears to strongly de- 
crease for T>3, it is rather constant for Sb and earlier 
galaxies. 

4. DISCUSSION 

4.1. The trend of central surface brightness with bulge 
luminosity 

Whether the central surface brightness of spheroidal 
components becomes brighter or fainter with increasing 
spheroid luminosity provides useful diagnostics on galaxy 
formation. The sign of this correlation results from the 
combined effects of dissipation (e.g., Kormendy 1989), 
feedback from star formation (e.g., Silk 2005), and dy- 
namical heating by merging black holes (e.g., F97). In 
mergers between unequal galaxies, it determines which 
component ends up populating the center of the merger 
remnant (Zurek et al. 1988; Balcells & Quinn 1990). Our 
trends of /io with luminosity (Fig.[5|) show that both the 
inwardly-extrapolated bulge /io, and /i measured at a 
fixed radius, become brighter with luminosity (i )3.1.3i 
Fig-Hfc, Fig.EJ. We have verified that fi measured at 
r = 20 pc from the center, and at a fixed fraction of 
i? c , specifically 0.2 R c , follow the same behavior. Fur- 
thermore, the slope defined in Figured by our bulges 
(brighter [j,q for more luminous objects) becomes steeper 
the closer to the center we measure fi. Jergen et al. (2000) 
have shown that, at the faint end of spheroid luminos- 
ity, the highly resolved, and also common, dwarf galaxies 
around the Milky Way and M31 have central profiles that 
do not show central excesses over their outer Sersic pro- 
file fits and have very low central surface brightness. The 
trend of brighter central surface brightness for higher lu- 
minosities is therefore well established for spheroids of 
vastly different luminosities, from giant ellitpicals to dEs 
(Fig.EK). 

Our findings contrast with the arguments of F97 who 
favor of a trend of fainter central surface brightness with 



increasing luminosity for ellipticals and bulges. F97 rea- 
soned that, if M32 resided in the Virgo Cluster, its 
exceedingly high central surface brightness would be 
smoothed by the PSF to values close to those of other 
low-luminosity Virgo E's. By the same argument, bulges 
and low-luminosity ellipticals at distances comparable to 
Virgo might have M32-like central profiles which, with 
sufficient spatial resolution, would yield /io values fol- 
lowing a faint-ward extrapolation of the trend defined 
by giant ellipticals, with Mb + 5 log/150 < —20.5 mag, 
in Figure [5Jd. 

Although Jerjen et al. (2000) have shown that the pic- 
ture described by F97 breaks down at the low-luminosity 
end, one could argue that the central galaxy surface den- 
sity, rather than bulge surface density, of some nucleated, 
intermediate-luminosity galaxies and bulges may follow 
the trend argued for by F97 (unresolved components [PS] 
reside in the centers of 58% of our bulges, see Paper III). 
Specifically, if the central surface brightness of the PSs, 
^o,ps 5 became sufficiently brighter for fainter PS (and 
hence bulge) luminosities, the trend between bulge lu- 
minosity and galaxy central surface brightness could in 
principle have an opposite slope to that shown in Fig- 
ure[5]for the trends between bulge luminosity and bulge 
central surface brightness. 

We find this unlikely. Both the central surface 
brightness of the underlying Sersic component (io,Scr ~ 
{Lk,Bui)°' 81 , eqn.O, and the PS luminosities (Lp$ ~ 
(£-/f,Bui) ' 91 , Paper III) rapidly become fainter toward 
lower luminosities, requiring exceedingly high densities 
for the unresolved components (over 5 magnitudes in sur- 
face brightness) in order to reverse the trend of Figured 
Boker et al. (2004), using HST/PC imaging, find that 
the sizes of nuclei do not correlate with nuclei luminosi- 
ties, hence the central surface brightness of the nuclear 
clusters Io.ps does not correlate with Lpg (for a sample 
of late- type spiral nuclei.) 

The compact elliptical galaxy M32 does provide a 
counter-example to the previous reasoning. The ground- 
based i?-band surface brightness profile of Kent (1987; 
see Graham 2002), and the color profiles from Peletier 
(1993), yield /i#(50pc)=14.64 mag. Given the absolute 
magnitude of M B = -15.74 for the M32 bulge 9 , M32 
lies within the sequence of bulges and low-luminosity el- 
lipticals in FigureO Yet, the HST profiles in F97 yield 
/z#(0.1")=8.52 mag, i.e., five magnitudes above bulges 
and ellipticals of the same luminosity. Note that the M32 
Sersic profile shape index derived by Graham (2002) is a 
modest n — 1.5, well within the range for bulges of sim- 
ilar luminosities; the high central brightness that makes 
M32 a prototype, compact elliptical galaxy really refers 
to a region of radius smaller than 37 parsec only (Gra- 
ham 2002 excluded radii smaller than 10" from his fit). 
This inner region is M32's 'unresolved source', using the 
term applied to our bulges. If the unresolved sources 
in bulges had similar structure to the M32 nucleus, they 
could reach very high central surface brightnesses despite 
their low luminosities. However, it is unclear whether 

9 We assume a bulge-to-disk ratio of \og{B / D) = 0.22 from 
Graham (2002); we use the apparent blue magnitude from LEDA, 
and a distance modulus of 24.43 (Jacoby et al. 1992). Ignoring the 
bulge-disk decomposition in determining the absolute magnitude 
would not affect our argument. 
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M32 is representative of spheroids such as those ana- 
lyzed in this paper; the bulge of M31 might be a more 
adequate model; it does not show such high central sur- 
face brightness. 

4.2. Bulge luminosity as a measure of the structure of 
disk galaxies 

Of the four sets of correlations involving bulge lumi- 
nosity, disk luminosity, bulge central velocity dispersion, 
and galaxy Hubble type T, relations involving bulge lu- 
minosity show the strongest level of significance. This 
was perhaps expected for bulge intrinsic parameters, but 
it is significant that disk parameters and bulge-to-disk 
ratios correlate with bulge luminosity as well. They also 
do so with the central velocity dispersion do, as a result 
of the correlation between oo and luminosity for bulges. 
In contrast, the dependencies with T are weak or non- 
existent. 

Such lack of scaling with T imply a stronger form of 
the statement from de Jong (1996) and Courteau et al. 
(1996) that the Hubble sequence is scale-free. These au- 
thors noted that bulge and disk scale lengths are corre- 
lated, and that their ratio R c /h is independent of Hubble 
type. We reproduce those results in our sample, as do 
HPG03 in theirs. Here, however, the lack of dependence 
on Hubble type is not restricted to length scales: the 
distributions of most bulge, disk, and bulge-to-disk pa- 
rameters versus T shows an absence of trends (Fig.0]) for 
our early-type disk galaxy sample. HPG03 find a similar 
absence of correlations for a later-type sample, selected 
with different criteria (galaxies belonging to the Perseus- 
Pisces supercluster) , and analyzed with a different (2D) 
code, suggesting that this behavior may be common for 
bulge-disk galaxies. 

Although the Hubble sequence is scale-free, galaxies 
themselves are not. The bulge luminosity (and velocity 
dispersion) correlate well with the other properties of the 
bulge, the disk, and the scaling of these two components. 

What structural parameters define the Hubble se- 
quence? Focussing on the early-to-intermediate disk 
galaxies covered by our sample (SO-Sbc), structural 
changes with T are only found in the latest, Sbc galax- 
ies (lower B/D, lower Iq, Bui /Iq, Disk an d marginally lower 
central velocity dispersions). For types SO to Sb, Hub- 
ble type must be dominantly given by the spiral pattern 
(Block & Puerari 1999; Seigar, Chomey, & James 2003), 
dust content, and star formation activity. Furthermore, 
the latter must not be affected by the spheroid luminosity 
or mass, nor by the relative sizes and brightness of disk 
and bulge. It appears that the luminosity and size trends 
commonly associated with the definition of the Hubble se- 
quence appear when extreme late-type galaxies are com- 
pared to early-types; but those trends are not gradual, and 
are absent among early-to-intermediate disk galaxies. 

That bulge luminosity strongly correlates with the 
properties of both the disk and the bulge-to-disk ra- 
tios may be expected from galaxy formation scenarios in 
which the bulge precedes the disk. This includes mono- 
lithic collapse models akin to the Eggen, Lynden-Bcll, 
&: Sandage (1962) model for the formation of the Milky 
Way, or CDM-based hierarchical galaxy formation sce- 
narios in which disks grow around already formed bulges 
of merger origin (e.g. Baugh et al. 1996; Abadi et al. 
2003; Sommer-Larsen, Gotz, & Portinari 2003). Our re- 



sults are probably also compatible with the inside-out 
galaxy formation model of van den Bosch (1998). The 
early formation of bulges does not necessarily imply a 
classical, slowly-rotating spheroid structure for bulges; 
if star formation timescales exceed the dynamical time, 
such bulges might show rapid rotation and some of the 
disky properties highlighted by Kormendy (1993). Are 
our results compatible with secular evolution models in 
which the bulge grows from instabilities in an already 
formed disk? The scaling of bulge and disk parameters 
is generally taken to support secular evolution models 
(Courteau et al. 1996; Zhang 2004). Here, we find that 
most bulge parameters show non-significant correlations 
with the disk i-T-band luminosity, suggesting that bulge 
properties (R c , mass, central density) do not know how 
massive is the host disk. Perhaps more puzzling, the disk 
parameters show more well-defined trends with bulge lu- 
minosity than with disk luminosity. Indeed, only one 
galaxy parameter, namely the disk scale length, corre- 
lates with the disk luminosity [Lk, Dis ~ h~ 2 (eqn. [8]). 
Explaining these relations, or lack of thereof, between 
galaxy structural parameters, will remain key challenges 
for secular evolution models. 

This difficulty is compounded by the evidence of stel- 
lar population ages. We show in Paper I that the bulge 
populations in the present sample are as old as cluster el- 
lipticals. Our measurements were performed away from 
the disk plane, hence ongoing or episodic star formation 
of secular origin could take place in inner regions of the 
disk. But old population ages are inferred for the nuclei 
of later-type spirals from NIR spectroscopy (Bendo & 
Joseph 2004). To be compatible with this evidence, sec- 
ular evolution must have been of little relevance for the 
growth of the bulge, or else confined to look-back times 
of order 10 Gyr for our galaxies: the situation essentially 
reduces to the van den Bosch model; and, timescales need 
to be known to verify that the term 'secular' applies. The 
situation may be different for bulges of later types, and 
of barred galaxies, which are not addressed in this pa- 
per (see, e.g., Fathi & Peletier 2003; Castro- Rodriguez & 
Garzon 2004); although, in late- type bulges, young pop- 
ulations are statistically associated with morphological 
disturbances, pointing to external, rather than internal- 
secular, processes (Kannappan et al. 2004). 

Hence, from the population information in Paper I 
and the structural information in this paper, it appears 
that bulges of unbarred, early- to intermediate-type disk 
galaxies cannot have significantly grown from disk in- 
stabilities in the last ~T0 Gyr. A small amount of star 
formation is probably occurring in almost every galaxy 
nucleus, as evidenced by dust in the very centers of bulges 
(e.g. Paper I), and by the presence of young stars (Free- 
man & Bland-Hawthorn 2002). Bar dynamics could be 
feeding such activity, but, for early- type disk galaxies, 
star formation level is too low to affect the general bulge 
population. Our conclusions on bulge formation/growth 
are consistent with the current understanding that many 
'mature' galaxies were largely in place by z = 1 (e.g. 
Brinchmann & Ellis 2000; Simard et al. 2002), and that 
early-type galaxies in the field became red between red- 
shifts z = 1 and 2 (Eliche-Moral et al. 2006). 

This research has made use of the NASA/IPAC Extra- 
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galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and 
Space Administration. This research has made use of 



the HyperLeda database. The United Kingdom Infrared 
Telescope is operated by the Joint Astronomy Centre on 
behalf of the U.K. Particle Physics and Astronomy Re- 
search Council. 
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